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Motivated by recent angle-resolved thermal conductivity experiment that shows twofold gap sym- 
metry in the high field and low temperature C phase in a heavy fermion superconductor UPt3, 
we group-theoretically identify the pairing functions as E\ u with the /-wave character for all three 
phases. Those are consistent with this observation in addition to variety of existing measurements. 
By using microscopic quasi-classical Eilenberger equation with this identified triplet pairing function 
under applied fields, we develop detailed studies of the vortex structures for three phases, includ- 
ing the vortex lattice symmetry, local density of states, and the internal field distribution. These 
quantities are directly measurable experimentally via SANS, STM/STS, and NMR, respectively. It 
is found that in the B phase of low H and low T the double-core vortex is stabilized over the sin- 
gular vortex. For the C phase thermal conductivity data are analyzed to confirm the gap structure 
proposed. We also give detailed comparisons of various proposed pair functions, concluding that 
the present Ei u with /-wave, which is analogue to the triplet planar state, is better than i?2u or 
Ei g scenarios. Finally we touch upon the surface topological aspects of Majorana modes associated 
with the E( u state of the planar like features. 



I. INTRODUCTION 

UPtji belongs to the first generation of the family of 
heavy fermion superconductors together with CeCu2Si2 
and UBei3 and has unique superconducting properties 
among them. Soon after the pioneering discovery^ of the 
double superconducting transition, it is founder— that the 
phase diagram in H vs T plane consists of three A, B and 
C phases. The A (C) phase is in high (low) temperature 
and low (high) field while the B phase is in the low T 
and low H . It is rather obvious to recognize that the 
order parameter (OP) must be multi-components. The 
main argument is centered on how to understand this 
phase diagram, or on what OP can describe it in a con- 
sistent manner— £. Now the splitting of the superconduct- 
ing transition temperatures T c \ = 550 mK and T c i = 500 
mK is generally understood due to a symmetry breaking 
field for otherwise doubly degenerate pairing state. The 
identification of this symmetry breaking field is still not 
settled yet, but it is considered as coming from the anti- 
ferromagnetic (AF) ordering at Tjy = 5 K^r— or coming 
from the crystal lattice symmetry lowering that occurs 
at a higher temperature^. 

The remaining question is to identify the OP symme- 
try. The central discussions are on what causes the OP 
degeneracy, either due to the orbital par t 6 ' 8 ' 13 ' 14 or the 
spin part of the OP— ~— . The former scenario has a fun- 
damental difficulty where the so-called gradient coupling 
term in the Ginzburg-Landau (GL) functional inevitably 
avoids the observed crossing of the two transition lines 
starting from T cl and T c2 , removing the C phase. There- 
fore the orbital scenario needs a fine tuning of the under- 



lying Fermi surface topology and the detailed structure 
of the orbital functioniili 6 .. Among various proposals the 
Em symmetry is regarded as a most possible candidate 
where d(k) oc z(k\ — fc 2 + 2ik x k y )k z . This state is time 
reversal symmetry broken and fourfold symmetric in the 
A and C phases. In the B phase, there exist a line node 
in the equator and two point nodes in the poles. 

On the other hand, the spin scenario ^ 12 ' 17- — is free 
from this difficulty, but a difficulty to qualitatively esti- 
mate the spin-orbit (SO) coupling remains because the 
spin scenario assumes the weak SO coupling in contrast 
with the strong SO coupling assumption adopted in the 
E^u scenario 6 ' 13 ' 14 ' 16 ' 22 ' 23 . This controversy is resolved 
experimentally because the Knight shif t 24 ' 25 starts de- 
creasing below T C 2 when H ~ 2 kG for H || c. This 
field H mt corresponding to the rotation of the d- vector- 
gives an estimate of the SO coupling strength in this sys- 
tem, justifying the classification scheme from the weak 
SO coupling, which never attains in the strong SO case 
where the d- vector is strongly tied to the underlying crys- 
talline axes via the orbital part in OP. 

The basic requirements of the possible pairing state 
realized in UPt3 can be summarized as follows: (1) The 
gap structure contains both horizontal line node(s) and 
point nodes as evidenced by power law behaviors in vari- 
ous directional dependent transport measurements, such 
as thermal conductivit y 26 ' , ultrasound attenuation ex- 



periment; 



,28-31 



and also bulk measurements, such as spe- 



cific heat 3 - ' 32 ' 33 , penetration depth^i, nuclear relaxation 
time^ 5 -, and magnetization^ 6 - experiments. (2) As men- 
tioned, the detailed Knight shift experiment^ 5 - shows the 
decrease of the magnetic susceptibility below T C 2, de- 



pending on the field direction and its strength. Then it 
is concluded that the d- vector contains the 6-component 
and c-component in the B phase for the hexagonal crys- 
tal. Upon increasing H (|j c) this c-component turns to 
be the a-component. (3) Finally according to the recent 
angle-resolved thermal conductivity measurement^ it is 
found that the twofold symmetric gap structure in the 
basal plane for the C phase, the full rotational symme- 
try in the B phase and the horizontal line nodes at the 
tropical position, not on the equator of the Fermi sphere. 

In view of those previous^ and recent experiments^ 
we come to a new stage to critically examine the pro- 
posed pairing states belonging to the orbital scenario; 
the singlet category E\^ and the triplet category 
t 13 ' 14 in addition to the so-called accidental degen- 
- and also belonging to the spin sce- 
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Previously we have tentatively chosen the orbital part 
from E\ u representation 1 ^— where the OP is non- 
unitary and also from E2 U representation^ in our spin 
scenario among the classified pairing functions in the 
abscnc o 41 i 42 and in the presence of the antiferromag- 
netism 4 ^. However, the precise orbital form is still to 
be determined. In view of this finding we are now in the 
position to identify the precise pair functions for all A, 
B, and C phases, based on the spin degeneracy scenario. 

The arrangement of the paper is as follows. We first 
classify the pairing function group-theoretically and ex- 
amine the possible symmetry realized in connection with 
various existing experimental data in Sec. II. The for- 
mulation is given in Sec. Ill based on the quasi-classical 
Eilenbcrgcr theory. With the identified pair function we 
investigate the vortex structures in the B phase, includ- 
ing the exotic vortex core structure, namely the so-called 
double-core structure associated with multi-component 
OP in Sec. IV. The present study here is a natural exten- 
sion of our GL framework for multiple OP— — to micro- 
scopic calculation based on Eilenberger theory, that pro- 
vides the quasiparticle (QP) structures in periodic vortex 
lattices. The vortex structure is also studied for the A 
and C phases in connection with the twofold gap func- 
tion in Sec. V. The final section is devoted to conclusion 
and future perspectives where detailed comparisons with 
other proposed pairing functions are made. We use the 
two notations (x, y, z) and (a, b, c) to denote the spatial 
coordinates interchangeably. 



II. CLASSIFICATION OF POSSIBLE PAIRING 
STATES 

In this section we first enumerate possible pairing 
states group-theoreticall y 41 i 42 ' 47 . The classified states 
are all so-called inert phase that is stable under a small 
change of the system parameters and differs from the 
previous classifications where we only treat the p-wave 
state a 41 i 42 . By the same procedure as in refs. |4l| andl42l. 
we can easily extend it to the /-wave pairing states with 
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FIG. 1. (Color online) Schematic phase diagrams under H || c 
(a) and H _L c (b). Orbital states are \ a {k) = ^21/%k a (5k 2 c - 
1) and \b(k) = ^/21/8kb(5k 2 — 1). (c) Gap functions of the 
three A, B, and C phases. 




FIG. 2. (Color online) DOS in the bulk for A and C phases 
|A(fc )| = A \k x (5k 2 - 1)| (solid line) and B phase |A(fc)| = 
Ao\/fci + ky\5ki — 1| (dashed line). 



the orbital angular momentum 1 = 3. 

Table I lists all possible /-wave inert states and its little 
group under hexagonal symmetry D e . Here the basis 
functions of the irreducible representations are also given 
in Table I. As mentioned before the required properties 
are all satisfied by the planar state f x lf lu + Tyl^ 1 " in 
E\ u . Namely in terms of the present context of UPt3, 
(ckh + bfc a )(5fc 2 — 1) where the unit vectors a, b and c in 
Dq denote the d- vector component. Note that this state 
is a natural extension of the planar state f x k x + T y k y 
in the p-wave state that is realized in thin films of the 
superfluid 3 He B-phase. 



A. Phase assignment and gap structures 

This Ei u state is assigned to each phase as the B phase 
(ckb + bk a )(5kc — 1), the A phase bk a (5k^. — 1), and the C 
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TABLE I. /-wave inert phase and its little group in Dg system. 
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A(e z ) = {u(e z ,e)9\0 <9< 2n}, C2 = {E+C 2z }, C' 21 = {E+C 21 }, D' 3 = {C 3 +C 21 C 3 }, 
m = {C 3 + C' 2 \ C 3 } , 11 D 2 = { 1 , C22 ^22 , C'21 u 2x , C' 2 \ u 2 y } , C 6 = {C 3 6 ( 2f ) , j = 0, 1 , ■ • • , 5} , 
cl = {Cj(^),j = 0,l,--- ,5}, nCl = {C J 6 u 2 6 j J = 0,l,--- ,5}, f M = t q M <r ;; (^ = x,y,z). 



phase ckb(5k^ — 1) as shown in Fig. [T] Thus the B phase 
is characterized by two horizontal line nodes at cos 2 9 = 
1/5, or 9 = 63.4° and 116.6° and two point nodes at the 
poles (see Fig.QJc)). The C (A) phase is the vertical line 
node at fcf, = (k a ~ 0) in addition to the two horizontal 
line nodes (see Fig.QJc)). Thus it explains naturally the 
twofold thermal conductivity oscillation when H rotates 
within the basal plane. Note also that this gap structure 
with two line nodes and two point nodes in the B phase is 
consistent with a variety of transport measurements and 
bulk measurements mentioned above. In Fig. [5] we show 
the density of states (DOS) N(E) for three phases. 

B. Unitary versus non-unitary 

It is seen from Table I that we list up the bipolar state 
given by f x lf lu + ifyl^ 1 ^ that is non-unitary and breaks 
the time reversal symmetry and that also can explain the 
phase diagram in H vs T. However, it is not consistent 
with the recent /xSR experiment*^ that negates the ear- 
lier claim that the time reversal symmetry is broken^. 
The gap structure in this bipolar state is characterized 
by | k x ± k y 1 1 5fc 2 — 1 1 in the B phase that has fourfold sym- 
metry. This contradicts the thermal conductivity exper- 
imental that indicated the rotational symmetry in the B 
phase. 



C. Other classified states and strong SO state 

The remaining states among the classified inert phase 
A 2 , B\, B 2 , Ei and E 2 are not accepted because of the 
following reasons: They do not provide the double tran- 
sition (A 2 , Bi, and B 2 ) 1 or do not explain the twofold 
symmetry in the C phase (i?2)- The other E\ u states in 
Table I are not appropriate except for the planar state, 
such as the polari, polar2, axial, pi, (3 2 , and 7 either be- 
cause they do not give the double transition (polari and 
polar2), or they fail to give an appropriate gap structure 
required (axial, /3i, p 2 , and 7). 

The E 2u c(fc 2 — k\ + 2ik a kt)k c classified in the strong 
SO case fails to explain the Knight shift experimen t 24 i 25 
and twofold symmetry in the C phased. Then we are 
left with only E\ u planar state with the /-wave character 
mentioned above. Note also that since E\ u with the p- 
wave character r x k x +if y k y has no line node, it has been 
excluded as a candidate from the outset. 



D. Dipole energy 

The dipole energy 

H D ^(3\k-d(k)\ 2 -\d(k)\ 2 ) k , (1) 
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depends on the combination of the spin and orbital 
states 5 ^, where (•••)* indicates the Fermi surface av- 
erage. The most favorable combinations by the dipolc 
energy in E\ u state are ak b (bk 2 — 1) and bk a (5k 2 — 1). 
In the C phase, the spin state b is selected by the AF 
ordering and accompanies the orbital state k a (bk 2 — 1) 
by the dipole interaction. In the B phase, the remain- 
ing orbital state k b (5k 2 — 1) has to appear with the spin 
state a to minimize the dipole energy. However, the pair- 
ing state in the B phase without the magnetic field is 



(bk a + ck b )(5kl - 1) actually (Fig. [TJ. Thus, the combi- 
nation of the spin and orbital states cannot be interpreted 
from only the dipole energy. 



III. QUASI-CLASSICAL EILENBERGER 
THEORY 



We start with the quasi-classical spinful Eilcnberger 
equation 5 -^—. The quasi-classical Green's function 
g(k,r,uj n ) is calculated by the Eilenberger equation 



-ihv(k) ■ Vg(fc,r,w„) 



{e/c)v{k) ■ A(r)} 
A(fc,r)t 



-A(fe,r) 
[ico n + (e/c)v(k) ■ A(r)} 1 



,g{k,r,u> Tl 



(2) 



where the ordinary hat indicates the 2x2 matrix in 
spin space and the wide hat indicates the 4x4 matrix in 
particle-hole and spin spaces. The quasi-classical Green's 
function is described in particle-hole space by 



g(k,r,w n ) 



g{k,r,w n ) if{k,r,u n )\ 
-if(k,r,u n ) -g{k,r,u n )J ' 



with the direction of the relative momentum of a Cooper 
pair k, the center-of-mass coordinate of the Cooper pair 
r, and the Matsubara frequency tu n = (2n+l)iik B T . The 
quasi-classical Green's function satisfies a normalization 
condition g 2 = —it 2 1. The Fermi velocity is assumed as 
v(k) = vpk on the three dimensional Fermi sphere. In 
the symmetric gauge the vector potential A(r) = (B x 
r)/2 + a(r), where B = (0, 0, B) is a uniform flux density 
and a{r) is related to the internal field B(r) = B + 
V x a(r). The unit cell of the vortex lattice is given 
by r = si(«i — 1*2) + S2W2 with —0.5 < Si < 0.5 (i = 
1,2), u\ = (a x ,0, 0), and M2 = (a x /2 7 a y} 0). In this 
coordinate, a hexagonal lattice is described by a y /a x = 
V3/2 or 1/(2^3). 

Spin triplet order parameter is defined by 



A(fc, r) = id(k, r) ■ crdy 



(4) 



with 



d(k, r) = aA a (r)(/> a {k) + bA b (r)4> b (k) + cA c (r)0 c (fe), 

(5) 

where & is the Pauli matrix. The self-consistent condi- 
tion for Ai(r) is given as 



A(fc,r) = N a Trk B T 

x ]T (V{k,k') [f(k',r,uj n ) + f(k',r,u n ) 

0<LU n <U c 



(0) 



I 

where No is the density of states in the normal state 
and uj c is a cutoff energy setting w c = 2QirkBT c with 
the transition temperature T c . We neglect the splitting 
of T c because that is appropriate for low temperatures 
even in the B phase. The pairing interaction V(k, k') = 
gcf>(k)(f)(k')* , where g is a coupling constant. The pairing 
functions 4>(k) and <fii(k) are chosen for each phase in 
UPt3 appropriately. In our calculation, we use a relation 



5^0 T c 



2nk B T 



E -• 



(7) 



0<OJ„ <LU C 



The vector potential for the internal magnetic field A(r) 
is also self-consistently determined by 



V x [V x A(r)] 

e 
c 



8TT-N 2TTk B T V (v(k) Im[g (k,r,io n )}) k , (8) 

n — ^ 

0<uj n <CJ C 



where g is a component of the quasi-classical Green's 
function g in spin space, namely, 
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.9o + 9z 9x - i9y 
fjx + ig y go - 9z 



We solve eq. J5]) and eqs. ^ and §5$ alternately, and 
obtain self-consistent solutions, under a given unit cell 
of the vortex lattice. The unit cell is divided to 41 x 41 
mesh points, where we obtain the quasi-classical Green's 
functions, Ai(r), and A(r). When we solve eq. © by 
the Riccati metho d 56 ' 57 , we estimate A;(r), and A(r) at 
arbitrary positions by the interpolation from their val- 
ues at the mesh points, and by the periodic boundary 
condition of the unit cell including the phase factor due 
to the magnetic field 5 ^— . In the numerical calculation, 
we use the units Rq ~ hvp / ' (2irk B T c ), B = ftc/(2|e|i?g), 
and Eq = irk B T c for length, magnetic field, and energy, 
respectively. By the dimensionless expression, eq. @ is 
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rewritten as where k = Bq / (EqV&xNo) = \/KW)7^ k G'l- We use a 

large GL parameter kql = 60 owing to UPt3. 

x [V x A(r)} 

Bo 

1 2T tt-^. By using the self-consistent solutions, free energy den- 

= — ^2 ^ ^ m r ' (9) sity is calculated by Luttinger-Ward thermodynamic po- 

c o<oj„<w c tentialSi as 



(5fi =A 



gN 



-TrA(fe,r)A t (fe,r) 



fc / r 

Tr 



V x A(r) - B 
Bo 



N J d\/irk B T (Rc[Tr{At(fc,r)(/ A (fc,r, W „)+£(fc,r, W „))} 



(10) 



0<Ul n <UJ 



where (• • • ) r indicates the spatial average. The auxiliary 
functions fx and / are obtained by the substitution of 

A A for A in eq. (|2|). This thermodynamic potential is 
relevant under large GL parameters and low magnetic 
fields because the replacement of the vector potential is 
not carried out. 

DOS for energy E is given by 

N(E) =(N(r,E)) r 

= (No (Re{go(k,r,u n )\iuj n -yE+i v ]) k ) r , (11) 

where r\ is a positive infinitesimal constant and N(r,E) 
is the local density of states (LDOS). Typically, we use 
r) = 0.0l7rfc B T c . For obtaining g Q (k, r, ui n )\ iUn ^,E+ir,, we 
solve eq. ([2]) with r\ — iE instead of uj n under the pair 
potential and vector potential obtained self-consistcntly. 

IV. B PHASE 

In the B phase, we take the pairing functions as <p = 
v / 2l78(fc a + h)(5k 2 c - 1), <f) b = ^2TJSk a {bk 2 c - 1), and 
<f) a = <f) c = y/21/8kb(5k 2 — 1), where one component of 
the d-vector is directed toward the &-axis and the other 
can rotate in the ac-plane. 

A. Double-core vortex lattice 

The pairing function in the B phase is similar to that 
in the superfluid 3 Hc B-phase, namely d(k) oc xk x + 
yk y + zkj£. Owing to the analogy with the 3 He B-phase, 
there is the possibility that the unconventional double- 
core vorte x 63 i 64 and v vorte x 65 ! 66 with a chiral core are 
stabilized against the conventional singular vortex. In 
fact, double-core vortex and v vortex are stabilized in 
the low and high pressure regions, respectively, in the 3 He 
B-phas o 63 ' 64 . Under our pairing function in the UPt3 B 
phase, the v (chiral core) vortex lattice is not stabilized 



self-consistcntly; however, there are two types of the self- 
consistent vortex lattice, namely, the hexagonal singular 
vortex lattice and the double-core vortex lattice. 

At T = 0.2T C under B = 0.055 to the c-axis, spa- 
tial variation of the pair potential amplitude for the 
singular vortex lattice and the double-core vortex lat- 
tice is shown in Figs. [3Ja) and EJb)-(d), respectively, 
where the total pair potential is defined by |A(r)| = 

J (Tr[A(fe, r)t A(fc, r)]/2)k- Since the pair potential is 
axisymmctric for the c-axis in the B phase, conven- 
tional singular vortices form a perfect hexagonal lattice 
(Fig. [3ja)). By contrast, a double-core vortex spon- 
taneously breaks the axisymmetry. Schematic struc- 
ture of the double-core vortex by cZ-vector is shown in 
Fig. EJc). The OP in the bulk is depicted by the blue 
(black) and red (gray) arrows which indicate d-vector 
with the orbital states A a = \j21/d>k a (hk 2 — 1) and 
Xb = \/21/8kb(5k 2 — 1), respectively. Along the 6-axis 
across the vortex center, the red (gray) arrow rotates in 
the ac-plane from the c direction far from the vortex to 
the — c direction at the opposite side of the vortex via 
the a direction at the vortex center. On the other hand, 
the blue (black) arrow directed toward the b direction 
is shortened as approaching the vortex center and finally 
vanishes at the vortex center; then, across the vortex cen- 
ter, the arrow is lengthened toward the — b direction to 
the initial length as far from the vortex. Thus, since d- 
vector can be modulated continuously across the vortex 
center, there is not the singularity where the total am- 
plitude vanishes. Instead of the singularity, the double- 
core with the small amplitude ~ O.6E0 exist as shown 
by contour lines in Fig.^b). This double-core has phase 
winding tt same as in the half-quantum vorte?* 6 ^. 

By the spontaneously broken axisymmetry, the double- 
core vortex lattice is distorted. The stable ratio between 
the height and base of the triangular lattice is a y /a x = 
V3/2.4, namely a base angle a = t&n~ 1 (2a y /a x ) w 55°. 
Each component of the double-core vortex lattice is also 
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shown in Figs. (3jc) and [3jd) for the amplitude of the 
bulk components |A&(r)| = |A c (r)| and that of the com- 
pensating component at the vortex cores |A a (r)|, respec- 
tively. The vortex core in the bulk component is slightly 
elliptic with the line of apsides along the a-axis and the 
compensating component is enlarged along the 6-axis. 
Since the vortex lattice tends to avoid the overlap of the 
vortex cores and that of the compensating component, 
the stable structure is fixed by the competition between 
them. 

At these low temperature and low magnetic field, the 
double-core vortex is more stable than the singular vor- 
tex. At high temperatures, however, the double-core 
vortex is unstable against the singular vortex because 
the compensating component tends to connect with the 
neighbor vortices along the o-axis by the extension of the 
coherence length. Similarly, under high magnetic fields, 
since the distance between the vortex centers becomes 
short, the double-core vortex is unstable. The Pauli- 
paramagnetic effect, which rotates the d- vector under 
H > ii/rot) is neglected in this calculation. Since the 
compensating component of the d-vector for b\ a + aXi, 
has to be directed toward the c-axis, the double-core vor- 
tex is unstable under high magnetic fields H > H lot also 
by the Pauli-paramagnctic effect. 

By the measurement of the small-angle neutron scat- 
tering (SANS), the perfect hexagonal lattice with a = 
60° is observed in the B phased. Thus, the observed vor- 
tices are conventional singular vortices. Since this exper- 
iment is carried out under the magnetic field H « H TOt , 
the double-core vortex may be unstable by the Pauli- 
paramagnetic effect. 



B. Local density of states 

There are clear differences in the LDOS between the 
double-core vortex lattice and the singular vortex lattice, 
which can be directly measured in scanning tunneling mi- 
croscopy/spcctroscopy (STM/STS). The LDOS for the 
double-core vortex lattice and the singular vortex lattice 
is shown in Figs. U and [5j respectively. The zero energy 
peak is expanded to the region between the double-core 
(Fig. 0Ja)) because the local OP aXb in this region has 
a line node in the ac-plane. Besides, there is the ellip- 
tic peak extending toward the a-axis in the LDOS at 
E = O.IEq (Fig. HKb)). By contrast, the singular vor- 
tex has an isotropic peak in the LDOS at E = and 
E = 0.1-Eo (Figs. 0(a) and E^b)). Spectral evolution 
of the LDOS near the vortex also differs between the 
double-core vortex lattice and the singular vortex lat- 
tice. Near the double-core vortex center, there is a sharp 
low energy peak especially along the a-axis (Figs. 0|c) 
and [He)) and along the 6- axis which is somewhat round 
(Figs. BJd) andQJf)). In the double-core case the OP 
two minima are situated just outside the center along 
the a-axis (see Fig.^b)), giving rise to sharp peaks at a 
finite energy as seen from Fig. B|e). On the other hand, 
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FIG. 3. (Color online) Spatial variation of the pair poten- 
tial at T = 0.2T C under B — 0.05Bo for the hexagonal sin- 
gular vortex lattice (a) and for the double-core vortex lat- 
tice with a w 55° (b)-(d). (a) Amplitude of the total pair 
potential |A(r)|, (b) amplitude of the total pair potential 
|A(r)| with the contour lines on 0.75i?o (solid lines) and 
0.85i?o (dashed lines), (c) amplitude of the bulk components 
|Ab(r)j = |A c (r)|, (d) amplitude of the compensating com- 
ponent at the vortex cores |A a (r)|, and (e) schematic spin 
structure of the double-core vortex around the core. 



the zero energy peak is varied to a bump away from the 
vortex core for the singular vortex (Figs. E](c)-[5jf)). It is 
noted that in the bulk region away from the vortex core 
the DOS oc \E\ 2 at low energy as seen from Figs. IHd) 
and EJd) because the gap structure is characterized by 
two point nodes and two line nodes. This is also seen 
from Fig. H 



C. NMR spectrum 

The double-core vortex lattice is also observed by the 
NMR measurement. In the NMR experiment, the reso- 
nance frequency spectrum of the nuclear spin resonance 
is determined by the internal magnetic field. The distri- 
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FIG. 4. (Color online) LDOS at T = 0.2T C under B = 0.05B 
for the double-core vortex lattice. Spatial variation of the 
LDOS at E = (a) and E = 0.1E (b). Spectral evolution of 
the LDOS from the vortex center (top) along the o-axis (c) 
and 6-axis (d). Details of the evolution near the vortex center 
are shown in (e) and (f), respectively. 
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FIG. 5. (Color online) LDOS at T = 0.2T C under B = 0.05B 
for the singular vortex lattice. Spatial variation of the LDOS 
at E = (a) and E = Q.lEo (b). Spectral evolution of the 
LDOS from the vortex core (top) along the a-axis (c) and 
b-axis (d). Details of the evolution near the vortex core are 
shown in (e) and (f), respectively. 



bution function is given by 

P(B) = J S[B - B c (r)]dr, (12) 

i.e., volume counting for B in a unit cell. This resonance 
line shape is called the "Redfield pattern" of the vor- 
tex lattice. In Fig. [51 we show the distribution functions 
P(B) for the singular vortex lattice (dashed line) and for 
the double-core vortex lattice (solid line). The distribu- 
tion function for the singular vortex lattice has a single 
peak at B = 0.049969i?o- This peak comes from the out- 



side of the vortex core shown by the contour lines in the 
left inset of Fig.[6j By contrast, the distribution function 
for the double-core vortex lattice has a double peak at 
B = 0.049962^0 and B = 0.049975B - The peaks at the 
low and high fields come from outside the vortex (solid 
line) and around the vortex (dashed line), respectively, 
shown by the contour lines in the right inset of Fig. [5] 
The distortion of the double-core vortex lattice gives a 
clear difference in the NMR spectrum. 
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FIG. 6. (Color online) Distribution functions of the internal 
magnetic field P(B) at T = 0.2T C under B = 0.05B for the 
singular vortex lattice (dashed line) and for the double-core 
vortex lattice (solid line). The height of P(B) is scaled so 
that J P(B)dB = 1. Inset: Spatial variation of the internal 
magnetic field B c (r) for the singular vortex lattice (left) and 
for the double-core vortex lattice (right) with the contour lines 
on the magnetic field situated at the peaks of P(B). 



V. C AND A PHASES 

In the C phase, we take the pairing functions as cf> = 
4> a = \/21/%kb(bk1 — 1) and <f)b — <f) c = 0, where the pair 
potential has one spin component. Note that the A phase 
is same as the C phase except that the a- and 6-axes are 
exchanged. 



A. Morphology of vortex lattice 

In this section, we discuss the deformation of the vortex 
lattice in the C phase under H \\ c to see the effects of 
the twofold gap function d(fe) cx ak^(hk\ — 1). Since the 
vortex cores are extended along the antinodal b direction, 
the height of the triangular lattice is enlarged along the 
b direction to avoid the overlap of the vortex cores. This 
variation of the vortex lattice is also same for the A phase 
by rotating it in the afr-plane. To find the stable vortex 
lattice, we compare the free energy among the triangular 
lattices with various ratios between the height and base, 
namely a y /a x . 

We show the stable vortex lattices under a low mag- 
netic field B = 0.02.Bo in Figs.[7(a)-[7(c) and a high mag- 
netic field B = 0.3B in Figs. [7(d)-[7(f). These figures 
arc also shown at different temperatures which are at a 
low temperature T = 0.2T C in Figs. [7(a) and [7(d), at 
an intermediate temperature T = 0.5T C in Figs. [7(b) 
and [7(e), and at a high temperature T = 0.7T C in 
Figs. [7(c) and [7^f ) . The triangular lattice is slightly 
distorted under the low magnetic field at the low tem- 
perature shown in Fig. [7(a). In this case, the ratio is 
a y /a x = 0.6-^/3, that is, a base angle of the isosceles tri- 



angular lattice is a = t&n^ 1 (2a y / a x ) ~ 64°. As tempera- 
ture increases, the stable structures of the vortex lattice 
are a y /a x = 0.65-\/3, namely a w 66° (Fig. EJb)) and 
a y /a x = 0.8^3, namely a 70° (Fig. EJc)). The vortex 
lattice is distorted drastically under the high magnetic 
field. At the low temperature, the stable structure is 
o. y /a x = 0.95Voj namely a w 73° shown in Fig. [7(d). 
Physically, the maximally distorted triangular lattice is 
a y /a x = y/Tf/2, namely a ~ 76° beyond which some of 
the nearest neighbors arc no longer nearest. Since the dis- 
tortion of the triangular lattice is near the limit even at 
the low temperature, the vortex lattice is hardly distorted 
by the increase of temperature as shown in Figs. [7(e) 
and [7(f) where a y /a x = V3; namely a ps 74°. Thus, 
the vortex lattice tends to be more distorted to avoid the 
overlap of the vortex cores under high magnetic fields 
and at high temperatures because the ratio of the radius 
of the vortex core proportional to the coherence length 
£ cx (T c — T) -1 / 2 to the distance between the vortices 
proportional to i? -1 / 2 increases. The deformation of the 
vortex lattice is summarized in Fig. [7(g). 

However, the regular hexagonal vortex lattice was ob- 
served in the A phase by the measurement of the SANS^. 
In this experiment, since the vortex lattice in the A phase 
is observed at the low temperature in which the B phase 
appears by the rapid cooling, the observed vortex lattice 
may change to the hexagonal singular vortex lattice in 
the B phase mentioned before. 



B. Local density of states 

The twofold symmetry of the gap function is revealed 
as shown in Fig [S] where the LDOS is displayed. The 
zero energy LDOS (Fig. [5(a)) is well connected between 
nearest-neighbor vortices along the a-axis which is the 
nodal direction. The LDOS at E = 0.LE (Fig. [5(b)) 
is also well connected between the nearest-neighbor vor- 
tices, and moreover, it has two peaks within a unit cell 
aligned along the 6-axis. Spectral evolution of the LDOS 
along the a-axis (Figs. [5(c) and [5(e)) forms the peak 
structure at a low energy near the vortex core. At few 
Rq ~ £ away from the vortex core, the spectra still in- 
crease from the zero energy so as to make a low energy 
peak. On the other hand, it has a rounded bump near 
the vortex core along the 6-axis (Figs. [5(d) and[5(f)). At 
few Rq away from the vortex core, the low energy spectra 
are almost flat. Far from the vortex core, their spectra 
become the same structure at the middle point between 
the vortex cores. It is noted that in the bulk region away 
from the vortex core the DOS cx \E\ at low energy as seen 
from Fig. [5(d) because the gap structure is characterized 
by one vertical and two horizontal line nodes. This is 
also seen from Fig. [5J 
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FIG. 7. (Color online) (a)-(f) Spatial variation of the stable 
pair potential amplitude |A a (r)|. (a) a ~ 64° under B — 
0.02B at T = 0.2T C , (b) a « 66° under B = 0.02B at 
T = 0.5T C , (c) a_« 70° under B = 0.02B at T = 0.7T C , (d) 
a w 73° under B = 0.3B at T = 0.2T C , (e) a w 74° under 
B = 0.3B at T = 0.5T C , and (f) a w 74° under B = 0.3B 
at T = 0.7T C . (g) Temperature dependence of the base angle 
a under B = 0.02Bo (solid circles) and B = 0.3Bo (open 
circles) . 



Field angle-resolved zero energy density of 
states 
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FIG. 8. (Color online) LDOS at T = 0.2T C under B = 0.02B 
in the C phase. Spatial variation of the LDOS at E = (a) 
and E = O.lBo (b). Spectral evolution of the LDOS from the 
vortex core (top) along the a-axis (c) and b-axis (d). Details 
of the evolution near the vortex core are shown in (e) and (f), 
respectively. 



We analyze the field angle-resolved thermal conductiv- 
ity experimental according to the identified pair function. 
In this experiment, since the temperature dependence of 
thermal conductivity obeys the Wiedemann-Franz law at 
low temperatures, the most significant effect on the ther- 
mal transport in the vortex state comes from the Dopplcr 
shift of the QP energy spectrum, E(p) — > E(p) — v s ■ p, 
in the circulating supercurrent flow v s . This effect be- 
comes important at such positions where the gap be- 
comes smaller than the Doppler shift term (|A| < v s -p). 



Since the magnitude of the Doppler shift strongly de- 
pends on the angle between the nodal direction and the 
magnetic field, the oscillation of zero energy DOS oc- 
curs. Consequently, thermal conductivity attains the 
maximum (minimum) value when the magnetic field is 
directed to the antinodal (nodal) directions^— . In this 
experiment, however, since heat current is injected along 
the c-axis, thermal conductivity cannot be compared to 
zero energy DOS directly. Then, we compare differences 



10 



\A a \/E D 



N/N 



0.2 0.4 0.6 0.8 1 1.2 05 1 1.5 2 




(c) 10 
5 

$ 



i i i i i 




i i i i i 


-10 -5 5 10 
a/R a 




-10 -5 5 10 
a/Bo 



C 

o 

> c 




-10 -5 5 10 -10 -5 5 10 

c/Ro c/Ra 



under various field directions, the field angle-resolved 
zero energy DOS is obtained shown in Fig. [TUtaL When 
field direction is rotated along the vertical line node from 
the c-axis (open circles), the zero energy DOS is dimin- 
ished because the low energy excitations from the tropical 
line nodes decrease. Within 45° < 9 < 60° and 120° < 
9 < 135°, the orientation of the vortex lattice changes, 
that is, the nearest-neighbor vortices are aligned along 
the 6-axis in 60° < 9 < 120° and next-nearest-neighbor 
vortices are aligned along the 6-axis in 9 < 45°, 135° < 9. 
By contrast, when magnetic field is rotated along the 
antinodal direction (solid circles), the zero energy DOS 
is almost constant because the QPs mainly come from 
the tropical line nodes under H \\ c and from the vertical 
line node under H \\ b. The difference of the zero energy 
DOS between the fields along the vertical line node and 
antinodc is maximum on the equator 9 = 90° because the 
horizontal line nodes are situated on the tropics. This 9- 
dependence of the difference of the zero energy DOS is 
consistent with the measurement of thermal conductiv- 
ity^, shown in Fig. fTUTbh Earlier spin singlet d-wave 
E\ g model (</> = vT5fct,fc c ) and spin triplet /-wave E^u 
model (cj) — <fi c = Vl05fc a fc(,fc c , <p a = <pb = 0)2, however, 
have two peaks of the difference of the zero energy DOS 
at 9 7^ 90° by the equatorial line node. 



VI. CONCLUSION AND PERSPECTIVES 



FIG. 9. (Color online) Spatial variation of the pair potential 
amplitude |A a (r)| and the zero energy LDOS N(r,E — 0) 
at T = 0.2T C under B = 0.05B . The magnetic fields are 
directed to the c-axis (a), b-axis (b), and a-axis (c). Left 
(middle) panels show the OP amplitude (zero energy LDOS). 
Right panels show a schematic view of line nodes seen from 
the field direction. 



of them when field directions are rotated along the ver- 
tical line node and antinode in the C phase. 

In this section, we assume a regular hexagonal vortex 
lattice a y /a x = or l/(2\/3). The stable orientation 

of the vortex lattice is determined by comparing the free 
energy calculated by eq. ([TP]). At T = 0.2T C under B = 
0.05i?o, the spatial variations of the pair potential ampli- 
tude and the zero energy LDOS are shown in left panels 
and middle panels of Figs. E]^a)-[9fc) , respectively. When 
the magnetic field is directed to the c-axis (Fig.[9£a)) or b- 
axis (Fig.lHJb)), elliptic vortex cores shrink to the vertical 
(H || c) or tropical (H || b) nodal directions. Under these 
magnetic fields, zero energy LDOS is mainly connected 
between nearest-neighbor vortices along the a-axis. On 
the other hand, under H || a (Fig.^c)), since the vertical 
node and the tropical nodes have a similar contribution 
to the QPs, vortex cores are hexagonal and zero energy 
LDOS is well connected among all nearest-neighbor vor- 
tices, giving rise to a rather round core profile in Fig.^c). 

By taking the spatial average of the zero energy LDOS 



In this paper we have classified possible pairing func- 
tions under the given crystalline symmetry D§h for the 
heavy fermion superconductor UPt3, which belong to the 
/-wave character with the triplet channel. Then we iden- 
tified a planar spin triplet state among them that maxi- 
mally fits to the existing experiments, in particular: 

(A) Various bulk thermodynamic measurements that 
indicate the line node(s) and point nodc(s) in the B 
phase. 

(B) The Knight shift experiment that shows the two 
field directions where the Knight shift decreases below 
T c , implying that the d- vector contains the 6-component 
and c-component at lower fields and the d- vector rotates 
from c to a at H TOt (|j c) ~ 2 kG. 

(C) The angle-resolved thermal conductivity measure- 
ment shows twofold gap structure in the C phase and 
rotational symmetry in the B phase. 

Those important experimental facts are all explained by 
the planar state with the /-wave channel, namely, (ckf, + 
bk a )(5k 2 c - 1). 

In order to check our proposed state, we made several 
predictions that are calculated by solving the microscopic 
Eilenberger equation with our planar state. The predic- 
tions include 

(1) The vortex structure in the C phase and likewise 
A phase exhibit strong distorted triangular lattice that 
varies as functions of T and H when H \\ c. This distor- 
tion is caused by the twofold gap structures in the A and 
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TABLE II. Candidate pair functions. 
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FIG. 10. (Color online) (a) Field angle-resolved zero energy 
DOS at T = 0.2T C under B = 0.05B - The field directions 
are rotated along the vertical line node (open circles) and 
antinode (solid circles) as schematically shown in the inset, 
(b) ^-dependence of the thermal conductivity normalized by 
the normal state value k„ (left axis) and the DOS difference 
normalized at 9 = 90° (right axis) along the vertical nodal and 
antinodal scanning for three possible gap functions in the C 
phase: 1. The present E\ u (kb(5k^. — 1)), 2. E\ g (kbk c ), 3. 
E 2u (k a kbk c ). The gap structures are sketched in the inset. 
The experimental data are taken from ref. l37l . 



C phases. The vortex morphology should be observed by 
SANS experiment. 

(2) While in the A and C phases the vortices are all 
singular, that is, the OP vanishes at the core because 
of the single component OP, in the B phase the vortex 
is non-singular characterized by the double core struc- 
ture. To check this complex vortex structure, we provide 
several signatures, such as the magnetic field distribu- 
tion probed as the resonance spectrum of NMR and the 
LDOS around a vortex core probed by STM/STS. 



A. Comparison with other proposed states 

There arc several proposals to identify the pairing sym- 
metries in UPt3. The E\ g and E^g are singlet and unable 
to explain the Knight shift experiment mentioned above. 
The E\ u with the p-wave character is not accepted be- 
cause there is no line node whose existence is firmly es- 
tablished through various thermodynamic experiments. 
The E2u which was regarded as most promising candi- 
date contradicts the observed twofold gap structure in 
the C phase. Table II summarizes the present status of 
various candidates. 



B. Remaining issues 

There remain several issues to be answered. 

(1) The d- vector rotation 

In order to understand the d- vector rotation phenomenon 
at H Iot ~ 2 kG for H || c, we need to take into account the 
magnetic field energy due to the anisotropic susceptibility 
in the superconducting state. In the absence of this effect 
the induced component a near the vortex core center in 
the double-core phase decreases as the vortex distance 
decreases. 

(2) The origin of the symmetry breaking 

To split T c into T c \ and T C 2, we need some symmetry 
breaking field. The good candidate is the AF order at 
TV = 5 K observed by neutron scattering. This is a high 
energy probe to catch the instantaneous correlation as a 
snap shot. Other low energy probes, such as NMR and 
/.tSR fail to observe the static AF order. Since a pre- 
cise correlation between the AF order and the T c split- 
ting under pressure, which simultaneously disappear at 
P ~ 3 kbar, is observed^, this is still a puzzling although 
we presented previously a scenario for this splitting due 
to AF fluctuations^. An alternative idea is to invoke 
the crystal symmetry lowering that is also reported bc- 
foreS^. 

(3) Pairing mechanism 

As mentioned in Sec. II the dipole energy alone cannot 
explain the stability of the pairing form (ckh + bk a )(5k^ — 
1). This special combination between the spin direction 
and the orbital form is rooted to the pairing mechanism. 
This is a main issue to be explore for future since now 
the pairing symmetry is determined in this paper. 

(4) Topological aspect 
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The identified pairing state (ckb + bk a )(5k^ — 1) is anal- 
ogous to the superfluid 3 He B-phase whose form is de- 
scribed as xk x + yk y + zk z realized in the bulk or the 
planar state xk x + yk y realized in thin films. Our state 
has the Majorana particles in the boundary as an An- 
dreev bound state albeit the line and point nodes exist 
in the bulk. This is interesting because recently Sato 
argues a possibility of the topological protection under 
the nodal gap 7 - 7 - The topological nature is discussed in a 
similar situation in connection with the superfluid 3 He 
A-phase where the Majorana particles exist in the point 
node gap 78 ' 79 . This topological aspect certainly deserves 
further investigation. Note that the double-core vortex 



does not contain the Majorana zero mode. 
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